Low-temperature phonon thermal conductivity of cuprate single crystals 
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The effect of sample size and surface roughness on the phonon thermal conductivity n p of 
Nd2CuC>4 single crystals was studied down to 50 mK. At 0.5 K, k p is proportional to y/A, where 
A is the cross-sectional area of the sample. This demonstrates that k v is dominated by boundary 
scattering below 0.5 K or so. However, the expected T 3 dependence of k v is not observed down to 
50 mK. Upon roughing the surfaces, the T 3 dependence is restored, showing that departures from 
T 3 are due to specular reflection of phonons off the mirror-like sample surfaces. We propose an 
empirical power law fit, to k p ~ T a (where a < 3) in cuprate single crystals. Using this method, we 
show that recent thermal conductivity studies of Zn doping in YBa2CusO y re-affirm the universal 
heat conductivity of d-wave quasiparticles at T — > 0. 

PACS numbers: 72.15.Eb, 74.72.-h 



To understand the pairing mechanism in a supercon- 
ductor, it is essential to know the symmetry of the or- 
der parameter. In this context, measurements of low- 
temperature thermal conductivity k, which probes the 
low-energy quasiparticle excitations, has emerged as a 
powerful probe of the order parameter in superconduc- 
tors. For conventional s-wave superconductors with a 
fully gapped excitation spectrum, the linear-temperature 
electronic contribution to thermal conductivity is zero at 
T — > 0, i.e. the residual linear term kq/T = 0. This 
can be seen in the single-gap s-wave superconductor Nb, 1 
multi-band s-wave superconductors MgB2 (Ref. 2) and 
NbSe2 r and two newly discovered superconductors C6 Yb 
(Ref. 4) and Cu^TiSea^ 

For unconventional superconductors with nodes in the 
gap, the nodal quasiparticles will contribute a finite 
Ko/T. For example, kq/T = 1.41 mW K~ 2 cm -1 in 
the overdoped cuprate T12201, a d-wave superconduc- 
tor with T c = 15 K£ and k /T = 17 mW K~ 2 cm" 1 for 
the ruthenate Sr2Ru04, a p-wave superconductor with 
T c = 1.5 Ki 7 - The fact that kq/T is universal^ in the 
sense that it is independent of scattering rate, allows for 
a measurement of the gap via^ 

k q k B 2 n v F 
1 6 c Ao 

This has been used to map out the gap as a function of 
doping in YBa 2 Cu 3 (y o and T1220l£ 

The measured thermal conductivity is the sum of two 
contributions, respectively from electrons and phonons, 
so that where n e /T is a constant as T —* 0. 

Therefore the key issue for these low-temperature ther- 
mal conductivity studies is how to extrapolate n/T to 
T = 0, i.e. to extract the residual linear term ko/T. This 
requires a good understanding of the phonon conductiv- 
ity K p (T). In the regime T — > 0, the phonon mean free 
path becomes limited only by the physical dimensions 
of the sample. At the surface of a crystal, the phonon 



may either be absorbed and reemitted with an energy 
distribution given by the local temperature (diffuse scat- 
tering) or it may be reflected elastically (specular reflec- 
tion). In the case of diffuse scattering, the phonon is 
reradiated in a random direction resulting in a temper- 
ature independent phonon mean free path Iq, given by 
the cross-sectional area A of the sample: Iq — 2\ZA/it. 
From simple kinetic theory, the conductivity of phonons 
is given by-ii 

K P = < v p > l T 3 , (2) 

where j3 is the coefficient of phonon specific heat, and < 
v p > is a suitable average of the acoustic sound velocities. 
The electronic linear term is then naturally extracted by 
plotting thermal conductivity data as n/T vs T 2 and 
interpreting the intercept as the residual linear term at 
T = 0, and the slope as the phonon contribution governed 
by Eq. (2). 

However, as the temperature of a crystal is reduced and 
the average phonon wavelength increases, a surface of a 
given roughness appears smoother, which may increase 
the occurrence of specular reflection and result in a mean 
free path which varies as some power of temperature, so 
that T 7 . We would thus expect a deviation from 

the diffuse scattering limit of T 3 temperature dependence 
for samples with sufficiently smooth surfaces. Such an 
effect has been previously observed in many studies of 
low-temperature phonon heat transport in high quality 
crystals, such as A1 2 3 ^ Siji 3 - KC1 and KBr^ii LiF^ 
and diamond^ 

For the high-T c cuprate superconductors, the extrap- 
olation of kq/T has been a controversial issue, partic- 
ularly in the underdoped regime where kq/T is small. 
Some authors fit their data to n/T = a + bT 2 below 
about 120 mK, assuming that boundary scattering of 
phonons only occurs below 120 mK and there is no specu- 
lar reflections i 17 i 18 Others consider that sample size and 
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FIG. 1: (color) In-plane thermal conductivity n of the un- 
doped insulating cuprate Nd2Cu04 below 0.5 K, for an as- 
grown single crystal before ("wide"; red squares) and after 
("narrow"; blue circles) it was cut along the length. By cut- 
ting, the width w was reduced by a factor of 5.6. At T — 0.5 
K, n is reduced in direct proportion to the reduction in \fw, 
as indicated by the solid line. 



FIG. 2: (color) In-plane thermal conductivity k of an as- 
grown single crystal of Nd2CuC>4 before ("as-grown"; full tri- 
angles) and after ("roughened"; empty triangles) its smooth 
mirror-like surfaces were roughened by sanding. As indicated 
by the horizontal line, roughening allows one to recover the 
diffusive regime at low temperature, characterized by a con- 
stant k/T 3 , at least down to 150 mK. 



specular reflections do affect n p , and fit their data to 
«/T = a + bT a below about 0.5 Disagree- 
ment on extrapolation procedure has fueled a debate on 
whether YBa2Cu30 y is a thermal metal below the dop- 
ing level of p = 0.05 , 17 ' 19 ' 20 and possible breakdown of 
the universal thermal conductivity in YBa2Cu307 and 
YBa 2 Cu 3 6 .5^ 

In this paper, we investigate the low-temperature 
phonon thermal conductivity of cuprate single crystals 
by measuring the insulating (undoped) parent compound 
Nd2CuC>4, for which Ko/T = 0. First, by reducing the 
sample width w, we show that k p is proportional to \fwi 
at 0.5 K, with t the unchanged sample thickness. This 
clearly demonstrates that k p is dominated by bound- 
ary scattering below 0.5 K. Secondly, we show a sam- 
ple with rough surfaces has a temperature dependence 
much closer to T 3 than the one with smooth surfaces, 
which is apparently caused by the reducing of specular 
reflections at the rough surfaces. Finally, based on this 
understanding of k p , we discuss how to extract kq/T in 
cuprate single crystals at finite doping (where kq/T > 0). 

Single crystals of Nd2CuC>4 were grown by a standard 
flux method. The as-grown samples are plate-like with 
very smooth surfaces in the a6-plane. In-plane thermal 
conductivity k was measured in a dilution refrigerator 
down to 50 mK using a standard one heater-two ther- 
mometer steady-state technique. Note that in zero field, 
magnons are gapped and the thermal conductivity of in- 
sulating Nd2CuC>4 only comes from phonons: k = k p 2^ 



Boundary scattering. — To study the boundary scat- 
tering of phonons, an as-grown Nd2CuC>4 single crystal 
with dimensions 1.5 x 0.90 x 0.086 mm 3 (length x width 
x thickness) was measured first. Then it was cut along 
the length so that its width w was reduced (by a factor 
of 5.6) to 0.16 mm, and measured again. Fig. 1 shows 
k vs T for the wide (before cutting) and narrow (after 
cutting) samples below 0.5 K. It is clear that the narrow 
sample has a much smaller k than the wide sample in 
this temperature range. Actually, at T = 0.5 K, n of 
the narrow sample is precisely \/5.6 = 2.37 times smaller 
than the wide sample. This means k is reduced in direct 
proportion to the reduction in \JA = y/wi, where t is 
the unchanged sample thickness. This is strong evidence 
that k is limited by sample size, and boundary scattering 
is unambiguously the dominant mechanism for phonon 
conductivity below 0.5 K. 

Specular reflection. - Having demonstrated that 
phonons are scattered by sample surfaces below 0.5 K, we 
proceed to study the effect of surface roughness on the 
temperature dependence of k. An as-grown Nd2CuC>4 
single crystal with dimensions 1.0 x 0.43 x 0.086 mm 3 
was measured first. Afterwards, both its mirror-like ab- 
plane surfaces were roughened by sanding, then it was 
re-measured. During sanding, the thickness of the sam- 
ple was reduced. The geometric factor for the roughened 
sample was set to be such that K(rough) = K(smooth) 
at high temperature (e.g. 50 K). In Fig. 2, the data of 
the smooth and roughened samples are plotted as k/T 3 
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FIG. 3: (color) In-plane thermal conductivity k of un- 
doped Nd2Cu04 (data from Ref. 21) and slightly overdoped 
YBa2Cu3C>7 (data from Ref. 18) single crystals, plotted as 
k/T vs T 2 below 175 mK. The dashed and solid lines repre- 
sent two extrapolation procedures: 1) k/T = a + bT 2 below 
120 mK and 2) k/T = a + bT a below 175 mK. The former 
yields a spurious residual linear term (kq/T > 0) in the in- 
sulator Nd2Cu04, and thus overestimates Ko/T in the super- 
conductor YBa2Cu3C>7. 
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vs T, to reveal the deviation from the T 3 dependence ex- 
pected if the phonon mean free path l p oc k/T 3 were con- 
stant, independent of T. l p is clearly not constant for the 
smooth sample, while for the roughened sample, rough- 
ening has made l p much more constant, at least down to 
150 mK. This is unambiguous proof that specular reflec- 
tion is important in these crystals. Below 150 mK, the 
phonon wavelength becomes long enough to average over 
the roughness and produce some specular reflections. 

Extracting Ko/T. — From Figures 1 and 2, it is clear 
that for typical cuprate single crystals with smooth sur- 
faces, phonon thermal conductivity will reach the bound- 
ary scattering regime below 0.5 K and k is expected to 
deviate from the standard T 3 dependence due to spec- 
ular reflections. Considering together previous studies 
of other conventional insulatorS ) 12 ' 13 i 14 ' 15 i 16 a reasonable 
way to extract kq/T of cuprate single crystals is to fit 
the thermal conductivity data to k/T = a + bT a . In Fig. 
3, we reproduce two sets of published data: one on un- 
doped Nd2Cu04 (from Ref. 21), the other on fully doped 
YBa 2 Cu 3 7 (From Ref. 18). For Nd 2 Cu0 4 , fitting k/T 
to a + bT a below 175 mK gives k /T = -0.004 ± 0.006 
mW K~ 2 cm -1 and a = 1.29 ± 0.04. In contrast, a fit to 
k/T = a + bT 2 below 120 mK gives k /T = 0.054 ± 0.003 
mW K~ 2 cm . For this insulator, the negligible kq/T 
obtained from the first fit is physically more reasonable 
than the finite kq/T from the second fit. This shows that 
a textbook fit overestimates Kq/T. 



FIG. 4: (color) Thermal conductivity of pure and Zn- 
substituted YBa2Cu30 H with (a) y = 7.0 and (b) y = 6.5 
(data from Ref. 18). The solid lines are fits to k/T = a + bT a 
below 175 mK, used to extract kq/T, seen to be universal in 
both cases (i.e. unchanged by Zn substitution). 



Applying the same two fits to YBa2Cu3C>7 (dashed and 
solid lines in Fig. 3) gives k /T = 0.162 ± 0.001 and 
0.115 ± 0.004 mW K~ 2 cm -1 , respectively. Based on 
our findings of Nd2CuC>4 (Fig. 3), the lower estimate 
(kq/T = 0.115 mW K~ 2 cm -1 ) is expected to be closer 
to the true value for YBa2Cu3C>7. 

Let us now turn to the controversial issue of whether 
thermal conductivity is universal in YBa2Cu30 y . In 
a d-wave superconductor, kq/T is due to nodal quasi- 
particles, and standard theory shows kq/T to be 
"universal" ; 22 i 23 i.e. independent of impurity concentra- 
tion, in the limit of weak scattering (F <C Aq). Upon 
increasing the scattering rate T so that it becomes a 
significant fraction of the gap maximum Aq, kq/T is 
expected to increase slightl y 22 i 24 (assuming the normal 
state is a metal). Experimentally, it was first observed 
in optimally-doped cuprates: YBa2Cu30j, (y = 6.9) as 
a function of Zn doping^ and Bi2Sr2CaCu20 , 8_ a ; as a 
function of radiation damage: 25 It was later verified in 
the layered nodal superconductor Sr 2 Ru04»i 

Repeating the original experiment, Sun et alr^ re- 
ported k data for a pure and a Zn-doped crystal of 
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YBa2Cu30j, with y = 7.0, reproduced in Fig. 4a. The 
main effect of Zn doping is to cause a dramatic suppres- 
sion of the slope of k/T as a function of T 2 . The authors 
extrapolate the data using a fit to k/T = a + bT 2 , re- 
stricted to a very small interval (70 to 120 mK). The 
value of kq/T thus extrapolated turns out to be 20% 
lower in the sample with 0.6% Zn. Based on this slight 
difference, Sun et al. go on to claim a breakdown of uni- 
versal transport in YBa 2 Cu30 y . 

However, the determination of k /T depends on the 
extrapolation procedure, as seen in Fig. 3. Using a fit 
to k/T = a + bT a below 170 mK for pure and 0.6% Zn 
doped YBa2Cu3C>7 (solid lines in Fig. 4a) yields Ko/T 
= 0.115 and 0.119 mW K~ 2 cm" 1 , respectively. This 
means that within error bars their data confirms the va- 
lidity of the standard theory of universal transport in this 
archetypal cuprate. In Fig. 4b, a similar result is found 
in underdoped YBa2Cu3C>6.5, where kq/T = 0.063 and 
0.058 mW K~ 2 citT 1 for pure and 0.6% Zn doped. 

Any variation in kq/T with Zn doping should be put in 
proper context: with 0.6% Zn doping, the inelastic scat- 
tering rate was estimated to increase by roughly a fac- 
tor 10 in YBa2Cu306.9^ Hence, the normal-state k/T in 
the zero-temperature limit should decrease by a factor of 
about 10. Therefore, the change in quasiparticle conduc- 
tivity measured in the superconducting state, if any, is 
seen to be at most a few percent of the change expected 
in the normal state. This is precisely what is meant by 
universal transport. 



In contrast to YBa2Cu30 a , the breakdown of stan- 
dard theory for La2_ x Sr 2 Cu04 (LSCO)^ appears to be 
fairly unambiguous (a rigid shift in k/T caused by a small 
change in scattering rate) and consistent with a prior re- 
port of breakdown in underdoped LSCOJ^ This is pre- 
sumably a consequence of the fact that suppressing the 
superconducting state in LSCO leads to an insulating 
(and magnetic) state, rather than the metallic state as- 
sumed by standard theory i 26 ' 27 

In summary, by studying the effect of sample size and 
surface roughness on the phonon thermal conductivity 
of Nd2CuC>4 single crystals, we show that the phonon 
heat conduction in cuprates is dominated by boundary 
scattering below 0.5 K. In as-grown (or polished) single 
crystals, specular reflection alters the T dependence away 
from the expected T 3 dependence. As a result, in no 
range of temperature down to T = is a T 3 fit to the 
phonon part of k appropriate. A better, but by no means 
exact, fit is obtained by allowing the power to adjust away 
from 3, towards 2, as found in conventional insulators. In 
the absence of a good theoretical treatment of specular 
reflection, this is probably the best one can do. 
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